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Abstract 

Background: The MADS-box transcription factors play fundamental roles in reproductive developmental control. 
Although the roles of many plant MADS-box proteins have been extensively studied, there are almost no functional 
studies of them in soybean, an important protein and oil crop in the world. In addition, the MADS-box protein 
orthologs may have species-specific functions. Controlling male fertility is an important goal in plant hybrid 
breeding but is difficult in some crops like soybean. The morphological structure of soybean flowers prevents the 
cross-pollination. Understanding the molecular mechanisms for floral development will aid in engineering new 
sterile materials that could be applied in hybrid breeding programs in soybean. 

Result: Through microarray analysis, a flower-enriched gene in soybean was selected and designated as 
GmMADS28. GmMADS28 belongs to AGL9/SEP subfamily of MADS-box proteins, localized in nucleus and showed 
specific expression patterns in floral meristems as well as stamen and petal primordia. Expression of GmMADS28 in 
the stamens and petals of a soybean mutant NJS-lOHfs whose stamens are converted into petals was higher than 
in those of wild-type plants. Constitutive expression of GmMADS28 in tobacco promoted early flowering and 
converted stamens and sepals to petals. Interestingly, transgenic plants increased the numbers of sepal, petal and 
stamen from five to six and exhibited male sterility due to the shortened and curly filaments and the failure of 
pollen release from the anthers. The ectopic expression of GmMADS28 was found to be sufficient to activate 
expression of tobacco homologs of SOCl, LEAFY, AGL8/FUL, and DEE. In addition, we observed the interactions of 
GmMADS28 with soybean homologs of SOCl, API, and AGL8/EUL proteins. 

Conclusion: In this study, we observed the roles of GmMADS28 in the regulation of floral organ number and petal 
identity. Compared to other plant AGL9/SEP proteins, GmMADS28 specifically regulates floral organ number, 
filament length and pollen release. The sterility caused by the ectopic expression of GmMADS28 offers a promising 
way to genetically produce new sterile material that could potentially be applied in the hybrid breeding of crops 
like soybean. 
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Background 

Many transcription factors (TFs) that control floral deve- 
lopment in model plants, such as Arabidopsis, have been 
identified and constitute a genetic "ABC model" for speci- 
f)^ing four-whorl floral organs: sepals, petals, stamens, and 
carpels [1,2]. Most TFs in the "ABC model" are so-called 
MADS-box transcription factors, which play fundamental 
roles in floral developmental control [3-5]. The ABC 
model proposes three classes of genes, termed A, B, and 
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C. A genes alone specif)^ sepals, A and B genes together 
specif)^ petals, and B gene expression in combination with 
a C gene specifies stamens. Carpel formation is dependent 
on C gene expression alone. In addition, the ABCD model 
was suggested when FLORAL BINDING PROTEIN II 
{FBPII), termed a D-class gene, was confirmed to deter- 
mine the ovule [6]. Although ABC genes are critical for 
floral organ identity, their expression levels are not suffi- 
cient to convert leaves into floral organs, indicating that 
other genes are required. Further studies in Arabidopsis 
have identified E genes {SEPALLATAI-4, previously 
named AGL2, AGL4, AGL9, and AGL3), which are 
required to specif)^ the identity of all four whorls of the 
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floral organs and floral meristem determinacy [7-11]. Ac- 
cordingly, the ABCDE model was proposed [12-14]. 

The SEPALLATA (SEP) or SEPAike genes have been 
described necessary, albeit redundantly, for the normal de- 
velopment of petals, stamens, carpels, and sepals [8,13,15]. 
It has been proposed that the SEP proteins could con- 
stitute higher-order complexes with A, B, or C proteins 
[8,13,15], and at least one SEP protein is present in more 
than half of the 106 multimeric complexes identified 
through a large-scale yeast three-hybrid screen [16]. SEP 
proteins are necessary for the formation of the transcrip- 
tion factor complexes that control ovule development 
[17]. A single sep mutant shows either a subtle or no 
phenotype, whereas the sepl/2/3/4 mutant displays inde- 
terminate flowers composed of only leaf-like organs; even 
sepal development is dependent on SEP function [10]. Al- 
though SEP4 is involved redundantly with the other SEP 
genes in the development of all floral organs, the sepl/2/ 
3/4 quadruple mutant exhibits the more extreme pheno- 
type of the loss of floral meristem identity compared to 
the sepl/2/3 triple mutants [10]. Additionally, the SEP3 
mutation alone leads to a partial transformation of the 
petals into sepals, suggesting that SEP3 plays more pre- 
dominant role in floral organ development than other SEP 
proteins [10]. SEP3 strongly activates the expression of the 
A gene API, B gene AP3, and C gene AG; these factors 
are the major protein interaction partners of SEP3, sug- 
gesting that SEP3 can activate the flower developmental 
program by enhancing the expression of its interaction 
partners [18]. The mutation of a rice 5'£P-like gene, 
OsMADS34/ PAP2, altered inflorescence morphology and 
increased the number of primary branches but decreased 
the number of secondary branches [19,20]. Moreover, 
osmads34 mutant displayed a decreased spikelet number 
and had lemma/leaf-like elongated sterile lemmas. Similar 
to the sepl/2/3/4 mutants of Arabidopsis, knock-down of 
the four rice 5'£P-like genes led to homeotic transfor- 
mation of all the floral organs, except the lemma, into 
leaf-like organs, suggesting the conservation of 5'£P-like 
genes in specifying floral determinacy and organ identity 
in both eudicots and monocots [21]. 

The phenotypes of transgenic plants constitutively ex- 
pressing SEP genes have been analyzed to further study 
their functions. Pelaz et al. [22] found that the overpro- 
duction of SEP3 in Arabidopsis led to early flowering 
and that the secondary shoots were transformed into 
solitary flowers. Furthermore, the ectopic expression of 
SEP3 was found to be sufficient to activate other B and 
C genes, and the constitutive expression of SEP3 and 
LEAFY, a floral meristem identity gene, produced ec- 
topic petals, carpels, and ovules outside of the floral 
context [23]. Zhao et al. [24] found that the ectopic ex- 
pression of a wheat SEPAike gene, TaMADSl, caused 
early flowering and altered the development of all floral 



organs, including sepals converted into leaf-like struc- 
tures and decreases in the petal and stamen number. 
The overexpression of SEP3 genes from tobacco and rice 
in Arabidopsis or tobacco promotes early flowering but 
does not significantly affect floral morphology [25,26]. 

Soybean {Glycine max [L.] Merr.) is an important glo- 
bal crop that provides protein and oil. Understanding 
of the processes that occur during soybean reproductive 
development and the identification of the genes respon- 
sible for developmental processes is important for soy- 
bean breeding. However, only a few reproductive-related 
genes have been identified in soybean, including GmAPl 
[27], GmNMH? [28], and GmGALl [29]. Additionally, 
none of these soybean MADS -box genes were shown to 
alter floral morphology. As an attempt to understand 
flower developmental process in soybean and genetically 
modify soybean floral morphology for breeding appli- 
cations, we identified 28 flower-enriched transcription 
factors in soybean through a microarray analysis [30]. 
Here, we report the functions of GmMADS28 which 
plays pivotal roles in the regulation of floral organ num- 
ber and petal identity. The sterility caused by the ectopic 
expression of GmMADS28 offers a promising approach 
to genetically produce new sterile material that could 
potentially be applied in crop hybrid breeding. Through 
the identification of more genes related with soybean 
flower, engineering soybean sterile materials applied in 
hybrid breeding programs will be possible. 

Results 

GmMADS28 is a member of Class E MADS-box genes 

Among 28 transcription factor genes predominantly 
expressed in soybean flowers [30], a soybean homolog 
(Gma.l7031.1.Al_at) of Arabidopsis AGL9ISEP3 was 
selected for further analysis in the present study. The 
microarray data were confirmed by real-time qPCR 
(Figure la). We named this gene GmMADS28 because it 
showed the highest similarity with Arabidopsis SEP3, First, 
we cloned the cDNA fragment of GmMADS28 via RT- 
PCR from soybean flowers. Through RACE-PCR, we as- 
sembled and cloned the full-length cDNA of GmMADS28, 
which is 1,026 bp in length and contains an ORE of 
732 bp. The sequence of GmMADS28 has been deposited 
in GenBank/EMBL under accession number AJ878424. 

Comparison of the GmMADS28 cDNA and soybean 
genomic DNA sequences suggested that GmMADS28 
contains eight exons and seven introns. By comparison of 
the exon-intron structures of GmMADS28 and other 
plant MADS-box genes, we found that GmMADS28 and 
other plant SEP genes share the same numbers of exons 
and introns; each exon size is highly conserved, whereas 
the intron size is divergent (Figure lb). The maize AG 
gene ZAG2 has seven exons, though the sizes of the first 
six exons covering the MADS- and K-do mains are similar 
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Figure 1 Characterization of GmlVlADS28. (a) Real-time qPCR analysis of GmMADS28 in soybean leaves, roots, flowers and pods, (b) Exon-intron 
structures of GmMADS28 and Lotus joponicos LjSEPS, Arabidopsis SEP3 and maize Z/\G2. The black and white boxes represent exons and introns, 
respectively, (c) Phylogenetic analysis of GmMADS28 and other plant SEP proteins. The neighbor-joining tree was constructed by MEGA 4.0 based on 
1,000 bootstrap replicates, (d) Subcellular localization of GmMADS28-GFP fusion protein. 



to those of GmMADS28, These results suggested that 
sequence alteration in the MADS -box genes during ge- 
nome duplications, particularly in the 3' region, might 
have resulted in divergent subfamilies of plant MADS -box 
proteins. 

The deduced GmMADS28 protein contains 243 amino 
acids and has a molecular mass of 27.9 kDa. Database 
searches suggested that GmMADS28 is homologous 
to Arabidopsis SEP3 (identity of 74% and similarity of 
83%), SEPl (identity of 59% and similarity of 71%), SEP2 
(identity of 60% and similarity of 71%), and SEP4 (iden- 
tity of 47%, and similarity of 64%). The alignment of 
GmMADS28 and the other plant MADS-box proteins 
showed that the amino acid sequences in the MADS 
domains were highly conserved. In addition to the 
MADS domain, the intergenic region (I region), K-box, 
and a C-terminal region (C region) are included in each 
MADS-box protein. Moreover, a potential phosphory- 
lation site, RQVTF, for calmodulin- dependent protein 



kinases and three common sequences, the "SEP motif I", 
"SEP/AGL6 motif, and "SEP motif IF [31,32], were found 
in all the analyzed SEP proteins, including GmMADS28. 

To address relationship of GmMADS28 to other plant 
SEP proteins, a neighbor- joining phylogenetic tree was 
constructed based on the alignment of amino acid se- 
quences from the plant SEP proteins (Figure Ic). The tree 
showed that GmMADS28 was grouped into the plant 
SEP3 subfamily. In particular, GmMADS28 was closer to 
Lotus japonicas LjSEP3 and Arabidopsis SEP3. These re- 
sults suggest that GmMADS28 is an E gene; based on the 
tree and highest sequence identity between GmMADS28 
and SEP3. 

GmMADS28 is localized in the nucleus 

To assess the subcellular localization of GmMADS28, we 
fused the full-length ORE of GmMADS28 to the green 
fluorescence protein (GFP) reporter gene under the con- 
trol of the CaMV35S promoter to generate the construct 
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35S:GmMADS28-GFP. The construct and the empty 
vector were transformed into onion epidermal cells via 
Agrobacterium-mediated transformation. Confocal im- 
aging of GFP fluorescence showed that the GmMADS28- 
GFP fusion protein was localized to the nucleus; in con- 
trast, the GFP fluorescence was distributed throughout 
the cells when the cells were transformed with GFP alone 
(Figure Id). 

Transcripts of GmMADS28 accumulate predominantly in 
reproductive organs, especially in petals 

The expression profiles of the GmMADS28 in various soy- 
bean tissues were investigated by semi- quantitative RT- 
PCR. Various soybean tissues, including leaf, flower, root, 
shoot apex, seed, stem, and pod, were collected. The ex- 
pression analysis showed that the GmMADS28 mRNA 
was clearly detected in the reproductive organs, including 
the flower, seed, and pod, but not in the leaf or root 
(Figure 2). Weak expression of GmMADS28 in the vegeta- 
tive shoot apex was also found (Figure 2). To further 
analyze the expression of GmMADS28 in flowers, we ex- 
amined the expression of GmMADS28 in the four whorl 
organs: sepal, petal, stamen, and carpel. As shown in 
Figure 2, GmMADS28 was clearly detected in all four or- 
gans and showed the highest expression in petals. 

Because GmMADS28 was highly expressed in deve- 
loping seeds (Figure 2), it is of interest to investigate 
whether GmMADS28 is involved in soybean seed develop- 
ment. Thus, we performed a semi-quantitative RT-PCR 
assay to analyze the time-course expression pattern of 
GmMADS28 during soybean seed development. The 
results showed that GmMADS28 was highly expressed at 
15 and 20 DAF (days after flowering) but was down- 
regulated thereafter, with the lowest level at 35 DAF; the 
expression of GmMADS28 increased gradually after 40 
DAF (Additional file 1: Figure SI). These results suggested 
that GmMADS28 may play some roles at the early and 
late stages of seed development. 

In situ localization of GmMADS28 transcript 

To further analyze the spatial distribution of GmMADS28, 
mRNA in situ hybridization was employed to detect floral 



mRNA expression. In the early stages of flower develop- 
ment, transcripts of GmMADS28 were localized in the 
floral meristem (Figure 3b), and in the stamen and petal 
primordia (Figure 3c-d). In the late stages of flower deve- 
lopment, transcripts of GmMADS28 were detected in the 
stamens, petals, and ovules (Figure 3e-f). 

GmMADS28 promotes early flowering 

To analyze the biological role of GmMADS28, we ectopi- 
cally expressed the GmMADS28 gene in tobacco (Figure 4a) 
and obtained 13 independent transgenic lines overex- 
pressing GmMADS28 confirmed by RT-PCR analysis 
(Additional file 2: Figure S2). We found that all positive 
35S:GmMADS28 plants flowered significantly earlier than 
the wild-type plants. As the plants have a fixed number of 
leaves when flowering under the same environment, a 
comparison of the flowering time was evaluated by the 
leaf number. The 35S:GmMADS28 plants produced 
an average of only 22 leaves when flowering, whereas 
the wild-type plants produced approximately 34 leaves 
(Figure 4b). The 35S:GmMADS28 plants also exhibited 
reduced plant height compared to the wild-type plants 
(Figure 4c). 

Ectopic expression of GmMADS28 affects morphology of 
floral organs 

The GmMADS28 over-expression lines exhibited re- 
markable morphological alterations in the floral organs. 
In wild-type plants, the number of petals, stamens, and 
sepals is five, but six 35S:GmMADS28 lines displayed six 
for each organ (Figure 5a-c, e-h). Interestingly, antisense 
expression of GmMADS28 in transgenic tobacco reduces 
the numbers of petals from five to four but the numbers 
of other three whorl organs are still five (Figure 5d). Be- 
sides, four 35S:GmMADS28 lines converted the stamens 
into petals and two converted the sepals into petals. The 
ectopic expression of GmMADS28 could convert the se- 
pals and stamens into petals (Figure 5i, j), implying the 
regulation by GmMADS28 in petal identity. 

Moreover, the sepals of the three 35S:GmMADS28 trans- 
genic lines were more similar to the carpels than to the 
sepals of the wild-type plant (Figure 6a-c). The inner and 



GmMADS28 



Ac tin 



Figure 2 GmMADS28 expression in different tissues of soybean by semi-quantitative RT-PCR analysis. Actin gene was used as the 
reference gene. 
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Figure 3 In situ hybridization of GmMADS28 in soybean developing flowers, (a-e) The longitudinal sections of flowers at differential 
developmental stages hybridized with an antisense (b-f) or sense (a) probe, (e) The cross section of flowers at the mature stage, fm, floral 
meristem; sp, sepal primordium; pp, petal primordium; stp, stamen primordium; cp, carpel primordium; st, stamen; pe, petal; ca, carpel; an, anther; 
o, ovule. Bars: 100 |jm. 



outer epidermal cells of the tobacco sepal, which are char- 
acterized by the presence of globular-tipped trichomes, 
stomata, and irregularly shaped cells, were used for scan- 
ning electron microscopy (SEM) analysis (Figure 6d, g). 
However, similar to the wild-type carpel, the inner and 
outer epidermal cells in the 35S:GmMADS28 sepals were 
regular rectangular-shaped cells (Figure 6e-f, h-i). These 
results suggest that GmMADS28 might involve in the 
carpel identity. 

35S:GmMADS28 plants are sterile due to physical 
alterations in floral structure 

The 35S:GmMADS28 plants were sterile. To investigate 
the factors causing sterility, we examined the morpho- 
logical structure of the floral organs. As shown in Figure 7, 
eleven 35S:GmMADS28 lines showed shortened and curly 
filaments, with the filament being shorter than the pistil 
(Figure 7c, d). The lengths of the filaments were measured 
before and after stretching. The filament length of 35S: 
GmMADS28 was significantly shorter than that of wild- 
type (Figure 7e), and the stamens could not reach the 
stigma due to such a shortened filament. The SEM ana- 
lysis on the outer epidermal cells of the filaments was per- 
formed to analyze the reason for the shortened filaments. 
Compared to the wild-type filaments, which contain nor- 
mal cells (Figure 7f), the 35S:GmMADS28 filaments only 
have cellular wrinkles, which may lead to curly filaments 
(Figure 7g). 



In addition to shorter filaments, most of the 35 S: 
GmMADS28 anthers cannot split normally, and the pol- 
lens cannot be released (Figure 8a-b). To reveal the forma- 
tion of the unopened anthers, the anthers from 35 S: 
GmMADS28 and wild-type plants were collected at four 
flower development stages, and the anther morphology of 
paraffin sections was compared using a microscope. It was 
found that there were no significant differences between 
35S:GmMADS28 and wild-type at the beginning stages of 
anther maturation (Figure 8c, stages 1 to 2). However, 
although pollenss were released with the maturation of the 
wild type (WT) anthers (Figure 8c, stages 3 to 4), the 35 S: 
GmMADS28 anthers could not split (Figure 8c, stages 3 
to 4) and the pollens were not released. 

GmMADS28 activates the expression of S0C1, LEAFY, 
AGL8/FUL and DEF 

Based on the function of GmMADS28 revealed by its con- 
stitutive expression, we studied the expression of 9 genes 
involved in flowering time or organ identity in the 35S: 
GmMADS28 leaves. It was observed that the tobacco 
homologs of the flowering time genes SOCl and LEAFY, 
A gene AGL8/FUL, and B gene DEF were more accumu- 
lated in 35S:GmMADS28 compared to wild-type leaves 
(Figure 9), suggesting that GmMADS28 might directly 
regulate the expression of these genes, thereby controlling 
flowering time and petal identity. 
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Figure 4 Ectopic expression of GmMADS28 promotes early 
flowering, (a) Selected transgenic plants at flowering stage, (b) The 
number of leaves of wild-type and 35S:GmMADS28 plants when 
flowering, (c) The plant height of wild-type and 35S:GmMADS28 
plants when flowering. The error bars represent SD (n = 10). 



GmMADS28 is involved in the conversion of stamens to 
petals in a soybean mutant 

The 35S:GmMADS28 transgenic tobacco plants exhi- 
bited the conversion of stamens to petals, which is also 
the phenotype of soybean mutant NJS-lOHfs [30,33] 
(Figure 10a). Compared to wild-type soybean stamens. 



the top areas of most stamens are converted to petal- 
like structure in NJS-lOHfs. Therefore, it is of interest 
to investigate the expression of GmMADS28 in the 
flower organs of this mutant. We found that expression 
of GmMADS28 in the stamens and petals in NJS-lOHfs 
was higher than in wild-type NJS-lOHff (Figure 10b, c). 
The expression of GmMADS28 was particularly higher 
in the mutant stamens (-50 folds) compared to wild- 
type (Figure 10c), suggesting that the up-regulation of 
GmMADS28 may play a critical role in the conversion 
of stamens to petals in this soybean mutant. 

GmMADS28 interacts with soybean homologs of S0C1, 
API, and AGL8/FUL 

As MADS-box proteins can form different complexes 
with different functions during floral development, we 
screened for interacting partners of GmMADS28 by yeast 
two hybrid assay using a cDNA library prepared from soy- 
bean flowers. A total of five proteins (GmSIPl ~ GmSIPS), 
including soybean homologs of AGL8/FUL, SOCl, and 
API proteins, an ENTH/VHS family protein [34], and 
an unknown protein, were shown to interact with 
GmMADS28 (Table 1, Additional file 3: Figure S3). These 
proteins may form complexes with GmMADS28 and 
function in the regulation of flowering time and floral 
development. 

Discussion 

Despite their similar sequences, plant SEP3 genes vary in 
their expression patterns. SEP3 is expressed throughout 
flower development, from the floral meristem to the fully 
developed floral organs [35]. The SEP3 genes are ex- 
pressed in the inner three floral whorls of all species, with 
the exception of Aranda AdOMl, which is not expressed 
in stamens and carpels [36]. GmMADS28 also showed a 
weak expression in the soybean shoot apices in addition to 
flowers and seeds, indicating that GmMADS28 might also 
participate in biological events in the shoot apices. In fact, 
Gerbera GhGRCD2 is also expressed in vegetative tissues 
[11]. Although GmMADS28 was expressed in all the 
reproductive organs analyzed, including the four whorl 
organs, it showed the highest expression level in petals, 
suggesting that GmMADS28 might play an important role 
in general flower development but is particularly involved 
in controlling the identity and development of soybean 
petals. 

To further explore the GmMADS28 function, we ecto- 
pically expressed the GmMADS28 gene in tobacco and 
found that GmMADS28 promotes early flowering. Many 
studies have shown that the overexpression of SEPS-hke 
gene leads to early flowering [23,25,37,38] . In Arabidopsis, 
the SEP3 protein is potentially incorporated into com- 
plexes with proteins involved in flowering time, such as 
SOCl, AGL24, SVP, and AGL15 [16]. The Arabidopsis 
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Figure 5 Phenotypes of 35S:GmMADS28 at the reproductive stage, (a, e) WT flower with five petals and stamens, (b, c, f) 35S:GmMADS28 
flowers with six petals and stamens, (d) 35S:antisense GmMADS28 flower with four petals, (g) WT with five sepals, (h) 35S:GmMADS28 with six 
sepals, (i) 35S:GmMADS28 converts sepal to petal, (j) 35S:GmMADS28 converts stamen to petal. The arrows in (f) indicate the stamens. Bars: 
1 mm in (e-f). 

I, J 



SOCl'like genes AGL42, AGL71, and AGL72 promote 
flowering in the shoot apical and axillary meristems [39]. 
Interestingly, our data showed that GmMADS28 acti- 
vates the expression of the tobacco homologs of LEAFY 
and SOCl, implying that the activated expression of 
these two genes may lead to early flowering. However, 
SOCl was shown to be directly down-regulated by SEP3 
in Arabidopsis [18]. These contrasting observations may 
be due to the different plant systems. Enhanced early flo- 
wering was also observed in Arabidopsis plants over- 
expressing both SEP3 and API genes compared to plants 
over-expressing API alone [22]. Similarly, we observed 
that the ectopically expression of soybean GmAPl also 
promotes early flowering in Arabidopsis [27] . These results 
suggest that GmMADS28 and API may be incorporated 
into the MADS -box protein complex to cooperatively 
regulate flowering time, which is also supported by the 



interaction of the GmMADS28 and API proteins by yeast 
two-hybrid assay. 

In Arabidopsis, SEP genes are involved in the develop- 
ment of all four-whorl floral organs and that con- 
stitutive expression of SEP3 promotes the homeotic 
transformation of sepals and inflorescence meristems 
into carpelloid structures, including ovules [23]. In the 
present study, a number of floral developmental alter- 
ations were observed in GmMADS28 transgenic plants. 
Interestingly, not all transgenic plants showed the 
same/similar phenotypes, which is likely because that 
GmMADS28 as an E-type gene plays a wide role in all 
most aspects of floral development and different expres- 
sion level and genome insertion location may influence 
the phenotypes. Even in one plant, not all the flowers 
showed the same phenotypes. For example, in the lines 
with increased number of floral organs, about half of 
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Figure 6 35S:GmlVlADS28 develops the carpel-like sepals, (a) WT sepal, (b) WT carpel, (c) 35S:GmMADS28 sepal. The inner (d, e, f) and outer 
(g, h, i) epidermal cells of W sepals and carpels and 35S:GmMADS28 carpel-like sepals were observed through scanning electron microscopy 
(SEM) analysis. Bars: (a-c) 1 mm; (d, f, g, \) 10 pm; (e) 5 pm; (h) 2 pm. 



the flowers altered organ number, while in the sterile 
lines, all the flowers are sterile. 

Ectopic expression of GmMADS28 promoted the deve- 
lopment of extra set of floral organs in transgenic tobacco, 
implying the critical and positive role of GmMADS28 in 
the formation of floral organs. Interestingly, the antisense 
expression of GmMADS28 decreased the petal number 
from five to four, which might have resulted from the de- 
crease of native SEP3 expression in the tobacco flowers. In 
contrast, Zhao et al. [24] found that the ectopic expression 
of a wheat SEP-like gene, TaMADSl, decreased the petal 
and stamen number, suggesting that different SEP proteins 
might play diverse roles in the regulation of floral develop- 
ment. Our data suggest that GmMADS28 is involved in 
the formation of four-whorl organs and is particularly re- 
quired for normal petal and stamen development. The A, 
B, C, and SEP proteins likely act as multimeric complexes 
to activate downstream genes [8,10,22,40]. The expression 
analysis of several ABC genes in 35:GmMADS28 leaves 



suggests that the ectopic expression of GmMADS28 itself 
could activate the expression of the DEF and AGL8/FUL 
genes. The multimeric complexes formed with activated 
ABC proteins and GmMADS28 may promote the produc- 
tion of extra organs in the transgenic plants. 

The ectopic expression of GmMADS28 mediated the 
conversion of stamens and sepals into petals in transgenic 
tobacco. In our previous study, we identified a soybean 
mutant [33] with a phenotype similar to 35S:GmMADS28, 
with the conversation of stamens into petals. Interestingly, 
the expression of GmMADS28 in the stamens of the soy- 
bean mutant was significantly increased compared to the 
wild-type plants, and the mRNA in situ localization ana- 
lysis proved the predominant expression of GmMADS28 
in the petal meristem of soybean. These results indicate 
that GmMADS28 plays a critical and positive role on the 
B function. Our data showed that ectopic expression of 
GmMADS28 activates the expression of B gene DEF [41] 
in tobacco leaves, suggesting that the ectopic expression 
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Figure 7 35S:GmMADS28 develops the shortened and curly filaments, (a, b) WT stamens can touch the stigma, (c, d) 35S:GmMADS28 
stamens are shortened and curly and cannot touch the stigma, (e) Comparison of filament length of WT and 35S:GmMADS28. The blue column 
indicates the length of 35S:GmMADS28 curly filaments while the red column the straightened filaments. The error bars represent SD (n = 25). 
(t g) The epidermal cells of WT and 35S:GmMADS28 filaments were analyzed by SEM. Bars: (a-d) 1 mm, (f) 3 \im and (g) 5 |jm. 
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35S:GmMADS28 



WT 1 




Figure 9 Expression of 9 MADS-box genes in 35:GmlVIADS28 
leaves. Semi-quantitative RT-PCR analysis sliowed tine expression of 
9 MADS-box genes involved in flowering time {SOCl, LFY) and floral 
organ development in tobacco leaves. These tobacco genes are the 
closest homologues of the corresponding Arabidopsis genes. The 
accession numbers for these genes are shown in Additional file 4: 
Table SI. 



of GmMADS28 can confer B function. Moreover, an A 
class gene, AGL8/FUL, was activated by GmMADS28. Re- 
cently, it was found that FUL-like genes also control flo- 
wering time and petal identity in opium poppy [42] . 

Carpel-like sepals were observed in some of the 35S: 
GmMADS28 plants and the constitutive expression of 
SEP3 resulted in the homeotic transformation of sepals 
and inflorescence meristems into carpelloid structures in 
Arabidopsis [23]. It has been described that both SEP3 
and AG are required for functional C activity [10,40] and 
that SEP3 could ectopically activate AG [23] in 35S:SEP3 
leaves. However, unlike Arabidopsis, the expression of AG 
was not significantly regulated by GmMADS28 in our 
study, suggesting that GmMADS28 might have an AG- 
independent role in the formation of carpel-like sepals. 

Controlling male fertility is an important goal for 
breeding hybrid crops, and male sterility is associated 
with the failure of pollination in addition to the lack of 
viable pollen or pollen activity. We observed that the 
35S:GmMADS28 plants are sterile, which was caused by 
two major factors, shortened and curly filaments and the 
failure of pollen release, which are new findings of SEP3 
function. At the beginning stages of anther maturation. 
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Figure 10 Expression of GmMADS28 in soybean mutant 
NJS-IOHfs. (a) The flower without petals of the mutant. The arrow 
indicates the conversion of the stamens into petals, (b) Real-time 
qPCR analysis of GmMADS28 expression in petals of the mutant 
NJS-IOHfs and those of WT. (c) Real-time qPCR analysis of 
GmMADS28 expression in stamens of the mutant NJS-IOHfs and 
those ofWT. The error bars represent SD based on three replicates. 



there were no significant differences between 35S: 
GmMADS28 and the wild type plants. The epidermal 
cells absorb water from the endothecium cells due to 
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Table 1 Soybean proteins interacting with GmMADS28 



Name 


Gene ID 


Annotation 


GmSIPl 


Glymallg36110 


ENTH/VHS family protein 


GmSIP2 


Glyma05g32820 


Unknown protein 


GmSIPS 


Glyma04g31800 


AGL8/FUL 


GmSIP4 


Glymal8g45780 


SOCl 


GmSIPS 


Glymal6gl3070 


APETALAl 



transpiration during anther maturation, which may pro- 
duce tension in the un-thickened, adjacent regions of 
the two pollen sacs. For normal anthers, further water 
loss during maturation causes the parenchyma tissue be- 
tween the pollen sacs to rupture easily under the tension 
developed by the endothecium, thus releasing the pollen. 
For 35S:GmMADS28 anthers, it is hypothesized that the 
endothecium cells are thickened and the shriveled con- 
nection of the cellulose can't pull out the easy-rupture 
region, the tension was not sufficient to split the anthers; 
therefore, the pollen grains can't be released. 

Unlike the role of GmMADS28, the overexpression of 
the SEP3 genes from Arabidopsis, tobacco, and rice in 
wild-type Arabidopsis or tobacco promotes early flower- 
ing but does not significantly affect floral morphology 
[22,25,26]. These results suggest that plant SEP3 pro- 
teins have common functions in the regulation of flow- 
ering time. In this study, we observed the roles of SEP3 
proteins in the regulation of floral organ number and 
petal identity. More importantly, the sterility caused by 
the ectopic expression of GmMADS28 offers a pro- 
mising way to genetically produce new sterile material 
that could potentially be applied in the hybrid breeding 
of crops like soybean. 

Conclusion 

The molecular regulations on reproductive development 
in soybean are largely unknown. Additionally, controlling 
the fertility is an important goal in plant hybrid breeding 
but is difficult in some crops including soybean. Through 
microarray analysis, a flower-enriched gene GmMADS28 
encoding a MADS -box transcription factor was cloned 
from soybean. GmMADS28 belongs to E-type gene and 
may play a wide role in reproductive development. It was 
observed that constitutive expression of GmMADS28 in 
tobacco caused a number of reproductive development 
changes, including early flowering, conversion of stamens 
and sepals to petals, increased numbers of sepal, petal and 
stamens and carpel-like sepals. In particular, ectopic ex- 
pression of GmMADS28 caused sterility due to the shor- 
tened and curly stalks and the failure of pollen release 
from the anthers. GmMADS28 is thus a potential target 
gene for engineering the male sterile plants. Moreover, 
GmMADS28 was found to activate the tobacco homologs 
of SOCl, LEAFY, AGL8/FUL and DEF and interact with 



soybean homologs of SOCl, API and AGL8/FUL pro- 
teins, which provides the new clues in understanding the 
functional mechanisms for E-type proteins in plant repro- 
ductive development. 

Methods 

Plant materials 

The soybean {Glycine max L. Merr. cv Jackson) seeds used 
in this study were obtained from National Center for 
Soybean Improvement, Nanjing Agricultural University, 
China, and germinated in the experimental field of 
Nanjing Agricultural University. Mature leaves (fially ex- 
panded), roots, flowers (including flower buds and mature 
flowers), and pods at 20 days after flowering (DAF) were 
collected for the real-time RT-PCR analysis. Mature leaves 
(fully expanded), roots, stems (including internodes and 
nodes), and shoot apices were collected separately at the 
stage of the sixth expanding true leaf for semi- quantitative 
RT-PCR. At the flowering stage, mixed flowers (including 
flower buds and mature flowers), sepals, petals, stamens, 
and carpels were collected. Seeds were collected at 15, 20, 
25, 35, 40, 45, and 50 days after flowering. 

Cloning of GmMADS28 

To amplify a fragment (Probe ID: Gma.l7031.1.Al_at) 
of GmMADS28, we first obtained the expressed se- 
quence tag (EST) sequence for Gma.l7031.1.Al_at by 
searching the Affymetrix website (https://www.affymetrix. 
com/) using the NetAffx tool. A pair of primers (for- 
ward, 5 '-GAGATGGGAAGGGGAAGAGT-3; reverse, 
5'-ACAAATTGGATATCATCCTG-3') was synthesized 
for RT-PCR and used with flower cDNA as the template. 
The PCR conditions were as follows: 0.5 \i\ flower cDNA 
was amplified in a 25 [A volume containing 2.5 [A PCR 
buffer with MgCb, 0.5 [A 20 mM dNTPs, 1 [A Tag poly- 
merase, and 0.5 [A each specific primer (25 mM). The 
PCR amplification was performed using a DNA amplifica- 
tion machine (MJ, USA) with an initial denaturation at 94° 
C for 5 min, 33 cycles of 94°C for 30 s, 55°C for 50 s, and 
72°C for 1 min, and a final 72°C for 10 min. The PCR 
products were gel-purified, cloned into the pGEM-T vec- 
tor (Promega, USA), and sequenced (Invitrogen, Shanghai, 
China). Rapid amplification of cDNAs (RACE) was 
employed with SMART RACE technology (Clontech, 
USA) to obtain the full-length cDNA for GmMADS28, To 
amplify the full-length cDNA of GmMADS28, the follo- 
wing pair of primers was used for RT-PCR: for- 
ward, 5'-GAGATGGGAAGGGGAAGAGT-3', and re- 
verse, 5 '-CAAGGAAGAGGCTAGCTAGG-3 '. 

Real-time RT-PCR 

Real-time qPCR was performed in an optical 96-well plate 
using a BIO-RAD iQ5 real-time PCR system (BIO-RAD, 
USA), as previously described [30]. The soybean Actin 
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gene (GenBank accession No. V00450) served as a refer- 
ence gene. The threshold cycle (Ct) values of the triplicate 
PCRs were averaged, and the relative quantification of the 
transcript levels was performed using the comparative Ct 
method. To verify the microarray data for GmMADS28 
expression, the relative quantification relates the PGR sig- 
nal of the target transcript in the roots, flower mixtures, 
or pods to that in the leaves. The fold change was deter- 
mined by the following formula: fold change = 2"^^^^ 

where AAGt = (Gt^^rget gene ~ ^^Actin) flowers, roots, or pods ~ 

(Gttarget gene " C^Actm)iea.ves' the analysis of GmMADS28 
expression in the soybean mutant, the relative quantifica- 
tion relates the PGR signal of the target transcript in the 
organs of the mutant to that of the wild-type plant. The 
fold change was determined by the following formula: 



fold change = 2" where AAGt = (Gt 



-target gene 



(Ct, 



-target gene ' 



Ct, 



Actin) wild - type plant* 



Semi-quantitative RT-PCR 

Semi-quantitative RT-PGR was performed as previously 
described [30]. As a control, a 683 bp PGR fragment of 
the constitutively expressed soybean Actin gene was 
amplified. The primers for GmMADS28 were as follows: 
sense, 5'-GAGATGGGAAGGGGAAGAGT-3', and re- 
verse, 5'-GAAGGAAGAGGGTAGGTAGG-3'. For the 
RT-PGR analysis of nine tobacco MADS -box genes, total 
RNA was prepared from the leaves of the wild type and 
five 35S:GmMADS28 plants. The tobacco EF-la gene 
was used as an internal control. The primers for each to- 
bacco MADS -box gene are shown in Additional file 4: 
Table SI. 



RNA in situ hybridization 

Longitudinal and cross sections of soybean flowers were 
prepared as previously described [27,43]. RNA antisense 
and sense probes were generated from a 281 bp frag- 
ment of the 3' region of the GmMADS28 cDNA labeled 
with digoxigenin. The RNA in situ hybridization was 
performed as previously described [27]. 

Subcellular localization 

To produce the GmMADS28-GFP construct, GmMADS28 
was fused in-frame to the 5 ' terminus of the green fluor- 
escent protein (GFP) reporter gene under the control of 
the GaMV35S promoter. The following primers was used 
for the GmMADS28-GFP construction: 5'-GAGATGG 
GAAGGGGAAGAGT-3' (forward), and 5'-GAAGGAA 
GAGGGTAGGTAGG-3 ' (reverse). The construct was 
then used in Agrobacterium-medrnted transient expression 
in onion epidermal cells. The cells were examined by 
confocal laser-scanning microscopy (Leica TGS SP2, 
Mannheim, Germany). 



Ectopic expression in tobacco 

To address the function of GmMADS28, the sense and 
antisense GmMADS28 sequence were cloned into the plant 
binary vector pBI121. The recombinant plasmids were used 
for tobacco {Nicotianatabacum cv. SamSun) transfor- 
mation via the leaf disk transformation method [44]. 

Cell morphological analysis 

Tobacco anthers were collected at flower lengths of 0.2, 
0.5, 1, and 3 cm, fixed with Garnoys Fluid (ethanol: glacial 
acetic acid = 3:1) for up to 24 hours at room temperature, 
and then preserved in 70% ethanol at 4°G. After embed- 
ding in paraffin, the blocks were trimmed as necessary 
and cut into 8-10 \\m sections using a LEIGA RM2135 
microtome. These sections were stained with safranin and 
stained again with fast green. After mounting the sections 
onto slides, the tissues were observed and photographed 
with a microscope (Leica DMLB). 

Yeast two-hybrid assay 

As the full-length GmMADS28 protein has transcriptional 
activation ability (data not shown), the cDNA sequence 
encoding GmMADS28AG lacking the activation domain 
(185-243) was used to construct the GAL4-BD fusion 
pDEST32-GmMAD528zlC as the bait plasmid and used to 
screen a yeast cDNA library prepared from soybean 
flowers. The positive clones were then verified by retrans- 
formation and subjected to |3-galactosidase assay. 

Availability of supporting data 

The sequence of GmMADS28 has been deposited in 
GenBank/EMBL under accession number AJ878424. 
The GenBank/EMBL accession numbers of MADS -box 
proteins used in the phylogenetic tree are: GRGDl 
(AJ400623), GRGD2 (AJ784156), SEP1/AGL2 (M55551), 
SEP4/AGL3 (P29383), SEP2/AGL4 (M55552), SEP3/AGL9 
(AF015552), GmSEPl (DQ159905), OsMADSl (L34271), 
OsMADS5 (U78890), ZMM6 (AJ430692), ZMM27 (AJ43 
0694), VvMADS4 (AF373603), S1SEP3 (BAD10945), LjSEP3 
(AY770397). 

Additional files 



Additional file 1: Figure 51. GmMAD528 expression during seed 
development. Actin gene was used as the reference gene. DAF: days after 
flowering. 

Additional file 2: Figure S2. The RT-PCR analysis of 35S:GmMADS28 
transgenic plants. P: Positive control, plasmid DNA; V\^: wild type plant; 
1-13: the 35S:GmMADS28 lines. 

Additional file 3: Figure S3. The yeast two hybrid assay identified five 
proteins (GmSIPl ~GmSIP5) interacting with GmMADS28. The cDNA 
sequence coding GmMADS28AC lacking the activation domain was 
cloned in pDEST32 and used as a bait to screen cDNA library prepared 
from soybean flowers. The transformation of pEXP™32/Krevl and 
pEXP™22/RalGDS-wt or pEXP™22/RalGDS-ml served as strong positive 
and weak positive controls while the transformations of pEXP™32/Krevl 
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and pEXP™22/RalGDS-m2 or pDEST32 and pDEST22 served as negative 
controls. 

Additional file 4: Table SI. The primer sequences used in gene 
expression analysis. 
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